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ABSTRACT 

The number of orphan radio afterglows associated with 7-ray bursts (GRBs) that should 
be detected by a flux limited radio survey, is calculated. It is shown that for jetted GRBs this 
number is smaller for smaller jet opening angle 0, contrary to naive expectation. For a beaming 
factor f^^ = (0^/2)~^ ~ 500, roughly the value inferred by Frail et al. (2001) from analysis of 
afterglow light curves, we predict that between several hundreds to several thousands orphan 
radio afterglows should be detectable (over all sky) above 1 mJy at GHz frequencies at any 
given time. This orphan population is dominated by sources lying at distances of a few hundred 
Mpc, and having an age of ~ 1 yr. 

A search for point-like radio transients with flux densities greater than 6 mJy was conducted 
using the FIRST and NVSS surveys, yielding a hst of 25 orphan candidates. We argue that 
most of the candidates are unlikely to be radio supernovae. However, the possibility that they 
are radio loud AGNs cannot be ruled out without further observations. Our analysis sets an 
upper limit for the all sky number of radio orphans, which corresponds to a lower limit f^^ > 10 
on the beaming factor. Rejection of all candidates found in our search would imply f^^^ > 100. 
This, and the fact that some candidates may indeed be radio afterglows, strongly motivate 
further observations of these transients. 



Introduction 



Our understanding of gamma-ray bursts (GRBs) has been revolutionized by the discovery of x-ray 



afterglows by BeppoSAX (e.g. [Costa et al. 1997| ), and the following detections of optical transients (e.g. 



van Paradijs et al. 1997). These efforts led eventually to a conflrmation of the cosmological nature of GRBs, 



through both direct redshift measurements (e.g. Metzger et al. 1997), and imaging of the host galaxies 
(e.g. Sahu et al. 1997). The enormous power released in the explosion implies extremely large compactness 
of the source and, therefore, most GRB models involve compact or collapsed objects (e.g. Goodman 
1986| , [Paczyhski 1986i [Eichler et al. 1989| , [Woosley 1993| , [Levinson fc Eichler 1993| [Paczyhski 1998| ; see 



Meszaros 1999 for a review). In spite of the impressive successes of expanding relativistic "fireball" models 



(e.g. Paczyhski fc Rhoads 1993, Katz 1994, Meszaros fc Rees 1997, Victri 1997, Waxman 1997a, 3ari, Piran 



fe Narayan 1998; see Meszaros 2002 for review), the precise nature of GRB progenitors remains unknown. 



The total energy emitted in the GRB explosion and the rate at which such explosions occur in the 
Universe are important keys to the understanding of the progenitors. The determination of these factors 
is complicated, however, by virtue of relativistic beaming. During the phase of 7-ray emission the fireball 
expands with a large Lorentz factor, P ^ 10^'^, so that a distant observer receives radiation from a conical 
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section of the fireball of opening angle ~ fO around the line of sight. Thus, estimates of total energy and 
rate based on 7-ray observations are highly uncertain. The first evidence that the fireball may be jetted was 



provided by radio observations of GRB970508 ( Waxman, Kulkarni fc Frail 1998 ). These observations imply 
a jet of relatively wide opening angle, which expands sideways and approaches sub-relativistic, spherical 
expansion « 1 yr following the GRB (Frail, Waxman & Kulkarni 200C). The analysis of radio observations 
during the sub-relativistic phase allowed to determine the total GRB energy to be « 10^^ erg. 

Jetted GRBs have been widely invoked to explain the optical light curves of the afterglow emission, 
most notably so for the source GRB990123 ( ^tanek ct al. 1999i [Harrison et al. 1999| ). Frail et al. (2001) 



have analyzed a sample of GRB afterglows with known redshifts. They find that most bursts are jetted, 



500. 



with jet opening angle 6 ^ 0.1 and average "beaming factor," defined as /j^ = {9 /2)^ , of 
They also find that the total 7-ray energy release, when corrected for beaming as inferred from the afterglow 
light curves, is narrowly clustered around 0.5 x 10^^ erg. This result, although somewhat model dependent, 
also suggests that the conversion efficiency of fireball kinetic energy to radiation is high, in agreement with 
the conclusion of Freedman & Waxman (2001), who derived total fireball energy using early X-ray afterglow 
data. The energy estimates of Frail et al. are in agreement with those derived by Freedman & Waxman, 



and also with the total energy derived for GRB970508 (Frail, Waxman & Kulkarni 2000) 



Using radio surveys to constrain GRB beaming, rate and energetics through the search for radio 
emission from GRBs that were not necessarily detected in 7-rays, has been proposed by Perna & Loeb 
(1998) , Woods & Loeb (1999) and by Paczyhski (2001). Perna & Loeb have suggested to search for 
such "orphan radio afterglow" emission from GRB jets pointing toward us, while Woods & Loeb have 
suggested to look for ~ lO'^ yr old, non-relativistic remnants. Paczyhski suggested, based on the radio 
observations of GRB970508, to search for emission from ^ 1 yr old nearby GRB remnants, which have 
undergone the transition from a relativistic jet to spherical sub-relativistic expansion, during which emission 



is isotropic. Given the hints for an association of GRBs with supernovae (Galama et al. 1998, Bloom et 



al. 1999 , Reichart 199£ ), it was suggested by Paczyhski to search for strong radio emitters among nearby 
supernovae. 



In this paper we use the model, proposed in (Frail, Waxman & Kulkarni 2000, Livio & Waxman 



200c ) for the long term radio afterglow emission of GRBs, to calculate the expected number of orphan 
afterglows in radio surveys. We find that the detectable orphan population is dominated by ~ 1 yr old GRB 
remnants which have just undergone the transition to sub-relativistic, spherical expansion. We determine 
the dependence of the expected number of sources on uncertain model parameters, derive a lower limit to 
the expected number of sources, and demonstrate that independent constraints on the beaming factor of 
early GRB emission can be imposed using radio surveys. We searched for point-like radio transients, by 
comparing the FIRST and NVSS wide field catalogues, taken a few years apart, looking for sources present 
in one catalogue and absent in the other. Our approach is different than that proposed by Paczyhski, 
namely we look directly for radio transients without relying on an association with nearby supernovae, 
since the link between GRBs and luminous supernovae is not well established: The only clear case, 1998bw, 
had an atypical (very weak) GRB, while there are three low redshift, z < 0.45, GRBs^ with no prominent 
supernova reported. 



In 



In § 2 we provide a brief description of the afterglow model and its main results used in later sections. 
3 we calculate the number of radio afterglows expected to be detected above a certain flux. In § 4 we 



see tttp:/ /www. aip.de^jcg/grbgen. html 
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present a search for point-like radio transients using the FIRST and NVSS catafogs, and discuss possible 
contamination by radio supernovae and AGN. The implications of our results are discussed in § 5. 



2. An outline of the afterglow model 



The basic model assumes an ultra-relativistic, jet like ejecta propagating into an ambient medium of 
density n — Iuq cm^^, slowing down and expanding sideways, ultimately becoming non-relativistic. As 
long as the jet Lorentz factor F is larger than the inverse of the jet's opening angle 9, it behaves as if it were 
a conical section of a spherical fireball. Once F drops below 0~^, the jet expands sideways, and its behavior 



deviates from that exhibited in the conical phase (Rhoads 1997, Rhoads 1999). After a transition stage, 
in which the jet expands sideways, the flow approaches spherical symmetry and becomes sub-relativistic 
(Waxman, Kulkarni & Frail 1998, Frail, Waxman fc Kulkarni 2000| ). The transition from a jet to a spherical 
sub-relativistic evolution takes place over a time scale, 

tsNT ^ 6 X 10\E5i/no)'^^ s (1) 

( Frail, Waxman fc Kulkarni 200C , Livio fc Waxman 2000| ). Here, £^51 is the total energy carried by the 
jet, in units of 10^^ erg, corrected for beaming, i.e. it is related to the isotropic equivalent energy _Eiso by 
E = (0^/2)£'iso = fbEiso, where 9 is the jet opening angle. The subscript SNT refers to the Sedov-von 
Neumann- Taylor self-similar solution that describes the sub-relativistic flow. 

After the transition, the power is dominated by synchrotron emission from a sub-relativistic fireball 
(e.g. Frail, Waxman & Kulkarni 2000). Denoting by (^s) the fraction of thermal energy behind the shock 
that is carried by electrons (magnetic field) , and assuming that the electrons are accelerated to a power-law 



energy distribution, drie/d'ye cx 7^ ^ with p = 2, it can be shown (Livio 
frequencies above the synchrotron peak frequency at t 



Waxman 2000) that the flux at 
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is given by (see also Appendix of Frail, Waxman & Kulkarni 2000) 
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Here, is the luminosity distance, i?o = c/Hq — IQ^^h ^ cm, and vq is the observed frequency measured in 



GHz. Analysis of GRB970508 radio data implies (Frail, Waxman fc Kulkarni 2000) no ~ 1, 1/3 



3. Expected number of orphan afterglow sources in radio surveys 

Consider a survey with a flux threshold denoted by f^mm, and denote by the time, in units of 
tsNxi^ + z), at which the radio flux of a GRB at a given luminosity distance, d^, drops below detection 
limit. From eq. (^ we then obtain 

r™ = 77(l + z)5/9(dL/i?o)-'°/^ (4) 



where 

ImJy J \0-3/ \0-3/ V^g 



h'"^' k) k) (-) «rxr. (5) 
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We shall assume, in what follows, that the total energy output and beaming factor of GRBs are 
standard and, moreover, that the fraction of total energy release that goes into driving the external blast 
wave generating the afterglow emission, is the same for all GRBs. We make this assumption in order to 
simplify our analysis, and since the distribution of the above parameters is not well known. It is important 
to emphasize at this point that, as mentioned above, the analysis presented in Frail et al. (2001) indicates 
that the distribution of GRB true (i.e. corrected for beaming) gamma-ray energies is narrowly clustered 
around E^^^i ~ 0.5, with a corresponding average beaming factor of < f^^ >=< {9'^/2)-^ >~ 500. These 
results also imply high, order unity, radiative efficiency, consistent with the result obtained by Freedman & 
Waxman (2001) using a different method, and implying that the true total fireball energy is a few times 
E^, i.e. narrowly peaked around E^i ~ 1. Thus, we expect corrections to our results, due to deviations of 
parameters from their standard values, to be small. Moreover, our results, given as explicit functions of the 
standard parameter values, can be easily generalized to the case where parameter values are drawn from a 
general distribution, by replacing these values with their corresponding averages. 

Let Ti be the time at which the radio source just becomes isotropic. Then for a given flux threshold, 
the maximum luminosity distance below which a substantial number of sources can be detected is obtained 
from eq. (|) for Tm = n: 

dLn.ax{Zmax) = i?0 (??M )'/'" (1 + Z^nax)'^^- (6) 

For a given choice of GRB parameters, a flux threshold fmin at the observed frequency vg, and value of h, 
the parameters rj and are fixed. Equation (^ can then be numerically solved for the assumed cosmology 
to yield the maximum redshift Zmax encompassed by the survey. 

We denote by fi the observed number of GRB events per year per Gpc'^ aX z = (which, assuming 



that the GRB rate follows the redshift evolution of the star formation, is estimated to be n ~ 0.5 (Schmidt 



2001)), and by <&(z) the redshift evolution factor. The expected number of radio afterglows between z and 



z + dz can then be expressed as. 



'^^"■^N,dl^{r^-nMz), (7) 



dz dz 

where dl = R^^cdt is the light distance in units of i?o, and 

A^o = f^H^RlhtsNT ^ 3.5 X lQ\bmfbr^h-^n{E^i/naf'\ (8) 

The total number expected, Nr, can be obtained by integrating eq (^ from z — to z — z^ax, with 
Zmaxiv/n) determined from eq. (§). 

As an example, consider a survey with a flux threshold of 5 mJy. Assuming ^ 3, h ~ 0.75 and the 
above choice for the remaining parameters, we find d^max — 0.2i?o = 0.8 Gpc, and a corresponding redshit 
Zmax — 0.2. At such Small redshifts cosmological effects can be neglected to a reasonable approximation, 
and the luminosity distance can be taken to be equal to the proper distance. In this case eq. (|^) can be 
integrated analytically. One then finds, using n = 10^^n_icm~^ 
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As seen, Nji depends only weakly on r^. Consequently, any uncertainty in the time at which the radio 
emission becomes roughly isotropic should not affect our result considerably. We have also given the result 
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in terms of the isotropic equivalent energy, Ei^o, for which we have chosen a normaHzation of 10^^ erg, since 
for the sample of GRBs with known redshifts given in Frail et al. 2001, one finds < -E^jso >= 3 x 10^^ erg 
and < e]%1 >6/ii= 5 x 10^3 erg. Since the true energy inferred from the observations, which is constrained 
by the isotropic equivalent energy Ej^iso measured in GRBs with known redshifts, is inversely proportional 
to the beaming factor, Nr oc . This means that the number of orphans expected to be detected in 
a flux limited survey decreases with increasing beaming factor, unless the survey is sufficiently deep to 
detect all sources. This is, of course, due to the fact that larger beaming factor implies, for given isotropic 
equivalent energy, a lower true energy, and hence a lower radio flux at the transition to spherical expansion, 
as demonstrated by Eq. (||). 

The number of orphan radio GRBs, Eq. (||), depends on the model parameters E, n, and 
S^B- Afterglow observations strongly suggest that is close to equipartition, i.e. ^ 1/3 (Freedman 
& Waxman 2001), and that E is narrowly peaked around £'51 = 1 (Frail et al. 2001). The main 
uncertainty is thus in the parameters n and ^s. The peak flux of a GRB afterglow, for an observer lying 
along the jet axis, is proportional to EisoVn-^B, and for typical luminosity distance of 3 x 10^* cm it 



W^/^Bno/lQ-^EisoM mJy ([Waxman 1997b|, pruzinov fc Waxman 1999], [Wijers fc Galama 1998"|) 



Observed afterglow fluxes generally imply ^btt-q > 10 ^ for i?iso,54 ^ 0.1, and values ^btiq ^ 1 arc obtained 
in several cases (e.g. [Waxman 1997b , Wijers fc Galama 199^ Panaitescu & Kumar 2001). Our parameter 



choice in Eq. (g), (^b/0.3) = 0.1 and no — 0.1, is therefore conservative. Even for such a conservative 
choice of parameters we expect more than a dozen radio afterglows to be detected by an all sky survey like 
the NVSS, provided the beaming factor is not much larger than that anticipated from other considerations. 
Note, that our choice h = 0.5 is also conservative: If the GRB rate does not evolve with redshift, rather 
than following the star formation rate, the z = rate is about 8 times larger, n « 4, and the expected 
number of radio GRB remnants is correspondingly larger by the same factor. Finally, the brightest source 
should have a flux of 

30(500/,)-^/^(n/0.5)^/^ (^||) (^1^^^^19/36^11/9^-1/2 (^^/3)-V30 ^jy. (iq) 

Since the flux in this phase declines roughly as {t/tsNT)~^ , such a source should remain above 1 mJy for at 
least several years. This motivates follow-up observations of the brightest transients found in a survey. 

The derivation of eq. (^) does not take into account the contribution from sources having a lifetime 
shorter than TitsNT- These sources are expected to be beamed, and so the fraction that can be observed 
(those directed along our sight line) should be of course correspondingly smaller. On the other hand, the 
radio flux emitted at times earlier than ntsNT is larger and, consequently, these beamed sources can be seen 
out to a larger distance. We may estimate the number of beamed radio sources anticipated to be detected 
by the survey using the following argument. The ratio of the number of sources having an opening angle 
~ 6 that are pointing in our direction, to the number of sources that had just become isotropic, is given by 
{2TTd'^ /4:Tr)[T{9)/tsNT] for a two sided jet, where T{0) is the observed expansion time to opening angle 9: 



T{6) cx 9"^ (e.g. [Rhoads 1999|) . Noting that the radio flux measured by an observer within the beam drops 



as t 1/3 during the expansion of the jet (e.g. Rhoads 1999 ), the flux of a beamed source is higher than that 



of a source that just became spherical (lying at the same distance) by a factor 9~^/^ . This implies that 
beamed sources can be seen out to a distance larger by a factor 6*^1/3 compared to those that just became 
spherical. Thus, the ratio of the number of beamed sources to those that just became spherical, observed 
above a certain flux, is {9~^^^Y {1:9"^ / '^'K)[r{9) /tsNr] oc 6*3. Thus, the number of detectable sources is 
dominated by sources which have already expanded to spherical symmetry. 



For deeper surveys cosmological effects cannot be ignored and must be taken into account. We have 
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solved eqs. (^) and (Q) numerically for different choices of cosmological parameters. For the range of 
parameters presently favorable we find that Nfi depends only weakly on the assumed cosmology. In the 
following we present the results obtained for h = 0.75, ^Im = 0.3, and J^a = 0.7. We also made the popular 
assumption that the comoving GRB density distribution traces the star formation rate, and invoked a 
strong redshift evolution for the later, as suggested by Madau & Pozetti (2000). Specifically, we take 
<I>(z) = (1 + z)^ in the following calculations. 

The solution of eq. (|^) is depicted in fig la, where the dependence of the maximum redshift Zmax on 
rj/Ti is plotted. Fig. lb shows the dependence of Nh/Nq on 77 for different values of r^. 




l0g(Tl/X,) 




Fig. 1. — (a) A plot of the maximum redshift below which sources can be detected above a certain flux, fumim 
versus r]/Ti. The parameters 77 and Ti are defined in the text. For fixed afterglow parameters rj cx jy^Jn ■ (b) 
The cumulative number of radio afterglows Nb. above a certain flux, normalized to A'o, versus ry, for different 
values of Ti. the normalization factor A^o is given by eq. in the text. The analytic solution obtained in 
the Euclidean case in eq. for t; = 3 is displayed (red line) for a comparison 

For the cosmology adopted above, the radio flux of a source at a redshift of z = 2 is given by 

and is smaller by a factor of 4 at 2 = 3. Thus for our above (conservative) choice of afterglow parameters 
a flux limit of < 1/j,J?/ is required to see all the orphans if the GRB population is dominated by those at 
z<2. 



4. A Search for NVSS/FIRST Radio Transients 
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4.1. Analysis 

We conducted a search for point-like radio transients in the FIRST radio catalog (White et al. 1997). 
We searched for compact FIRST radio sources that do not appear in the NVSS radio catalog (Condon et al. 
1998). The sample we obtain this way is incomplete, since it rejects afterglow sources that were detected 
by the NVSS, and were bright enough to remain above 6 mJy during the FIRST observation. We examine 
this incompleteness factor further below and show that it is smaller than 5 percent. As the FIRST and 
NVSS surveys were obtained in different VLA configurations, the comparison was done carefully. We use 
this comparison to put an upper limit on the number of GRB radio afterglows. 

The FIRST is a 20cm radio survey conducted using the NRAO VLA in B-configuration, with 5" 
beam-size. The final maps have a typical noise RMS of 0.15 mJy. This value allows 5cr detection of 1 mJy 
sources. The survey's resolution enables an astrometric accuracy better than 1" (90% CL). The FIRST 
radio survey has begun in 1994, and when completed it will cover 10,000 deg^, mostly in the north Galactic 
cap. In the search described below, we used the FIRST 2001, Oct 15 version, covering 8,565 deg^. This 
catalog has 771, 076 sources, hence, a density of about 90 sources per deg^. 

The NVSS 20cm radio survey (Condon et al. 1998) was conducted between 1993 to 1996 and covers 
all the sky north of declination —40 deg. It utilized the VLA D- and DnC-configurations, with a circular 
beam-size of 45". This is significantly larger than the median angular size 10") of faint extragalactic 
sources. With this beam-size, the faintest NVSS sources have an astrometric accuracy better than 7" (68% 
CL). The images have RMS brightness fluctuations of ^ 0.45 mJy beam"^, that allows for 5a detection of 
^ 2.5 mJy sources. The NVSS 99% completeness level is about 3.4 mJy. With density of p ^ 60 sources 
per deg~^, it has total of ~ 2 x 10^ radio sources. 

Radio sources could be intrinsically variable. Therefore, using the NVSS 99% completeness of 3.4 mJy 
as our threshold, could introduce a large number of variable sources (e.g., that were detected in the FIRST 
survey, but were just under detection in the NVSS) into our transients list. In order to compare the two 
catalogs, we need to choose a flux limit that rejects most of the variable objects. We plot in Figure |^ all 
FIRST point-like sources (above 4 mJy) with NVSS counterparts. The solid line shows the location of 
sources for which the FIRST flux equals the NVSS flux, and the horizontal dashed line shows the NVSS, 
99% completeness flux level (i.e., 3.4 mJy). Based on this plot we chose a threshold of 6 mJy, for which the 
scatter around the solid line does not produce large number of false transients sources. 

In the following we assumed that all the FIRST images, were taken after the corresponding NVSS 
images. This assumption will be examined more carefully later. Unless a GRB afterglow happens in a radio 
galaxy, it is expected to: (i) be a point source; and (ii) be absent from the NVSS. Since radio galaxies are a 
tiny fraction of the overall galaxy population, these assumptions hold for most of the afterglow population. 

We selected all the FIRST point-like sources with peak flux above 6 mJy (21.5%). From this subsample, 
we selected all the unresolved objects (4.3% of the 6 mJy subsample). However, some of the resolved FIRST 
objects could actually be point sources, classified as resolved due to measurements errors. Adopting the 
FIRST size error criterion (White et al. 1997), for 9.2% of the FIRST sources above 6 mJy it is not possible 
to reject the hypothesis that they are point-like at the 99.7% CL. Therefore, the number of sources above 
6 mJy for which the point-source hypothesis could not be rejected at the 3(t level is 2.1 {— 9.2/4.3) times 
larger than the number of point sources we used in our sample. We ended with a catalog of 7181 FIRST 
point-like sources with peak flux above 6 mJy. For these sources, we searched the NVSS for counterpart 
within 30" of the FIRST position. Although the NVSS astrometry is better than 7", the use of such a high 



- 8 - 




Fig. 2.— FIRST vs. NVSS peak flux for all FIRST sources with NVSS counterparts. The solid hue shows 
the location where the FIRST/NVSS flux ratio equals 1, the horizontal dashed line shows the NVSS 99% 
completeness flux level of 3.4 mJy, and the vertical dashed line shows the 6 mJy cut. 

threshold was necessary because in cases of two (or more) nearby point-like radio sources, the NVSS could 
detect only one elongated source, with its center offset by more than 7" from one of the FIRST point source 
positions. The search yields 110 point-like FIRST sources with peak flux above 6 mJy for which there is no 
NVSS counterpart. 

We examined the POSS-II E images and the FIRST and NVSS radio maps of all these transient 
candidates. We found that: 2% result from holes in the NVSS maps; 14% are artifacts due to un-cleaned 
aliases in the radio maps; 61% are multiple sources for which the NVSS found an elongated source centered 
more than 30" from one of the FIRST radio positions; and for the rest (23%, 25 sources) we could not 
reject the hypothesis that these are radio transients or variable radio sources of some kind (e.g., AGNs, 
SNe, GRBs, Pulsars). A large fraction of these 'sources' are possibly: radio artifacts; NVSS non-detections 
due to locally high backgromid (e.g., in the vicinity of bright radio source); or physical phenomena other 
than GRBs (e.g., AGNs, SNe). If we change the flux limit from 6 to 10 mJy, the number of candidates 
decreases from 21 to 6. 

The following table lists the 25 candidates along with their basic properties. 

For each of these objects, we conducted a search for known counterparts within a radius of 5' in the 
NED database and within 1' in the HEASARC database. Our findings are listed in the table. 

2575 deg^ of the FIRST survey were taken before the completion of the NVSS survey. Therefore, we 
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reject all the sources in this area (i.e., 22.2° < 5 < 42.5°). Omitting this strip (PSR J1652+2651 is in this 
strip), the area covered by our search is therefore 5990 deg^, and the number of sources is: 16 and 3, for 
the 6 mjy and 10 mJy flux thresholds respectively. The upper 95% Poisson CL on these numbers are: 26 
and 9 respectively. We consider these numbers to be conservative upper-limits on the number of possible 
GRB afterglows. 

Finally, we have to correct our upper limits for the various completeness factors involved: (i) Missing 
point sources due to FIRST measurements errors. As we showed above, a very conservative estimate for 
this number is Cmissing = 9.2/4.33 = 2.1; (ii) Area completeness (e.g., the fraction of the sky covered by 
our survey) Carea = 41252/5990 = 6.8868; (iii) Given The crude NVSS resolution, the probability that 
a transient source will be detected up to r = 90" (twice the NVSS FWHM), from a known un-related 
NVSS source, and therefore could be confused with it, is P(< r) = 1 — exp {—npr^) = 11%. Therefore, 
the source confusion completeness factor is: Cconf = 1-11; (iv) NVSS completeness - For about 2% of the 
FIRST sources there are no corresponding NVSS images (holes in the catalog), hence we set Ccomp — 1-02. 
(v) Incompleteness due to our search criteria - afterglow sources detected by the NVSS that were bright 
enough to remain above 6 mJy during the FIRST observation are rejected by our search. The number 
of such sources depends on the average time separation between the NVSS and FIRST observations. 
Denoting by t^ep the time interval between two consecutive observations of a typical afterglow source, 
eq. (^) implies that the flux during the flrst observation was larger than during the second observation 
by a factor of {tsep/tsNT + 1)^^^", assuming the source was near maximum flux during the first detection 
(i.e., its age was roughly Isnt)- Consequently, the fraction of bright afterglow sources in the NVSS 
catalog that may be detected above 6 mJy by the FIRST at time tgep after 1996, must be smaller than 
5{tsep) = {tsep/tsNT + The Smallest time separation between the catalogs in the strips we consider 

is about tmin = 6 months, and the largest is about tmax = 4.5 years. Taking the FIRST detections to be 
uniformly distributed in time, integrating the fraction S{tsep) from imin to tmax yields a total fraction of 
less than 5% that our search rejected. Thus, we set Csearch = 1.05. Taking together all these factors, we 
can put a conservative 95% CL upper limit of 26 x Carea x Cmissing x Cconf X Ccomp X Csearch = 447 for the 
all-sky number of GRB afterglows with peak specific-flux above 6 mJy (or 154 for the 10 mJy threshold). 



4.2. Radio Super novae 



The characteristic behavior of radio supernovae (RSNe) (see [Weiler et al. 2002| for a recent review) 



appears to be similar to the expected behavior of radio afterglows (RGRBs). In particular, the radio light 



curve of type II RSNe rises on a characteristic time scale of ~ 100 days (Weiler et al. 1998), comparable to 
tsNT (see eq. [|l|), and has an overall similar shape to that of RGRBs following the peak. This is not at all 
surprising since both are produced by the same mechanism, namely a blast wave, albeit under somewhat 
different conditions. This similarity would render any attempt to distinguish RGRBs from RSNe difficult 
and, therefore, motivates a careful examination of potential differences between these two classes of objects. 

We identify two features that may allow one to distinguish RGRBs from RSNe. Firstly, the radio flux 
is anticipated to decline smoothly with a t^^ dependence during the sub-relativistic phase in RGRBs, while 
the radio light curves seen in RSNe are often much more complex. However, some RSNe do show a smooth 
decline, and those may be confused with RGRBs. Secondly, the absolute peak radio luminosity of RSNe 
may differ considerably from the radio luminosity of the brightest RGRBs. For a typical GRB the energy 
injected into the blast wave is E ~ 10^^ erg, and with the above choice of parameters: n — O.lcm^'^, t; = 3, 
S_B = 0.03 and = 0.3, the RGRB luminosity density at 5 GHz is 7 x lO^^erg s^^Hz-^ at a time of sa 0.5 yr 
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after the explosion. This luminosity is ~ 100 times larger than the typical peak luminosity (of detected) 



RSNe, of order lO^^erg s ^Hz ^ ( Weiler ct al. 1998 ), and several times brighter than the brightest known 
RSNe, SN1988Z, which had a peak flux density of 2 x lO^^erg s'^Hz-^ This suggests that in a flux limited 
survey RSNe and RGRBs should have a distinct redshift distribution, with RSNe having preferentially 
smaller redshifts. 

We estimate below the number of RSNe expected to be detected in a search of the type described in 
this paper for RGRBs, and compare their expected redshift distribution to that calculated for RGRBs. 
In order to make an accurate prediction for RSNe, their (peak) luminosity function and the shape of the 
corresponding light curves must be known. Unfortunately, due to the small sample of RSNe presently 
available (12 type II and 4 type I), and the lack of a supernova model that is capable of predicting the 
distribution of RSN properties, the determination of the peak luminosity function is highly uncertain. We 
therefore provide below a crude estimate based on the Weiler et al. (1998) compilation of RSNe radio data. 
Based on Figure 3 of Weiler et al., we adopt a RSNe peak luminosity function with equal number of RSNe 
per logarithmic peak luminosity interval over the 5 GHz peak luminosity range Li = lO^^erg s-^Hz-i to 
L2 — lO^^erg s~^Hz~^. In addition, we invoke a correlation between 5 GHz peak luminosity Lp and peak 
time Tp of the form, L.p = 102*(Tp/103day)i-4erg/sHz, as demonstrated in Weiler et al. (1998). These 
assumptions hold, based on the available data, for type II RSNe. The properties of Type I RSNe are less 
certain. However, since type II RSNe dominate the (detected) RSNe population, and since the type I RSNe 
peak luminosities are close to the typical type II peak luminosities, we expect our estimate below of the 
type II RSNe "background" rate to provide a (crude) estimate for the total RSNe background. 

Assuming that the RSNe flux drops as Tp/t for t > Tp, the number of 5 GHz sources on the sky above 
a certain flux density /i/^min at any given time can be readily computed: 

^^^^^ = 93h!^7ZO^^^^^"^^^'=^'- ^^^^ 

Here, fifjSNe is the RSNe rate per unit volume, T2 is the peak time corresponding to peak luminosity L2, 
and RfisNe is the maximum distance at which a RSN can be observed, 

/ L2 \'/' /L2/102WsHz\'/' , , 

RRSNe ^ =41 ^ j^L^ Mpc. (13) 

Since the typical distances of detectable RSNe are ~ 10 Mpc, we have assumed Euclidean geometry for the 
calculation above. It is also straight forward to calculate the average distance of detected RSNe, 

93 /La/lO^Serg/sHzy/^ 



< R >RSNe^ ^^RSNe = 13 ^ , "\ MpC. (14) 

304 V /2^,min/5mJy / 

These distances should be compared with those obtained for RGRBs: for the conservative GRB parameters 
quoted above, the maximum and average distances of RGRBs detectable above 5 mJy at 5 GHz are 
Rgrb = 140 Mpc and < R >grb~ Rgrb/3 = 46 Mpc respectively. 

Using the above equations, and taking fijfSNe = 1-2 x 10^ Gpc^^yr^^, corresponding to the total rate 



of core-collapse SNe ( Filippenko 2001 ), we find Nj^sNe = 3(/iy,min/5mJy)~'^/'^. Comparing this result with 
the estimates for RGRBs given in §2, we find that the number of RGRBs detected in a survey of the type 
described here should be at least comparable, and most likely significantly larger, than the number of 
detected RSNe. The characteristic distances of detectable RGRBs should be > 4 times larger than those 
expected for RSNe. 
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The anticipated proximity of RSNe in the FIRST survey impHes that the host galaxy should be easily 
detectable. In fact, even very faint galaxies should be detectable above the POSS-II limit out to a distance 
of 100 Mpc. This means that all the radio transients listed in table 1 for which no optical counterpart is 
detected are unlikely to be RSNe, as even dwarf hosts should be clearly evident. 



4.3. Radio loud AGNs 

Radio loud AGNs are known to exhibit large amplitude variations over a broad range of timescale 
and frequencies, with flare durations ranging from minutes to years. The characteristic durations of radio 



outbursts (with the exception of ID variability) are between weeks to years (Valtaoja et al. 1999, Vcnturi 



et al. 2001), with amplitudes that occasionally exceed a factor of 10. Thus, it is conceivable that some of 
the transients in our selected sample are radio loud AGNs. Indeed, our initial list of 110 transients contains 
a number of such sources. 

Most of the candidates in table 1 have no optical counterpart brighter than Mji = 21. What is the 
likeliness that some or most of them are in fact AGNs? One way to characterize the radio loudness is 
in terms of the radio-to-optical luminosity ratio, R, which is commonly defined using radio and optical 
luminosities at some fiducial rest-frame wavelengths. Existing samples of AGNs exhibit a distribution in 
radio loudness that depends on selection criteria. The distribution of R in the FBQS ( White et al. 200C| ) 



for example peaks at i? ^ 10, with very few sources (0.3 %) having R > 10"^. The distribution of radio 
selected quasars has a median at R ^ 10'^ (Wadadekar & Kembhavi, 1999). It is not clear at present 



whether this is a result of a selection effect, and whether there exists a population of extremely radio loud 
AGNs. However, based on our current knowledge, we would naively expect that deep enough observations 
at optical wavelengths would be able to ultimately reveal the underline AGN. future X-ray detections of 
some of our candidates will also imply that these are likely AGNs. A precise determination of the minimum 
value of R for the objects in table 1 requires a knowledge of their redshifts, which are not available. If R 
is defined as in Wadadekar & Kembhavi (1999) then a rough estimate yields R > 300 for those objects. In 
addition, they should have no extended structure above ~ ImJy on scales larger than a few arcseconds. 

Another potential difficulty for future surveys may arise due to the presence of faint, BL Lac type 
sources. If a population of radio BL Lacs exist for which the continuum flux associated with the nuclear 
activity drops below that of the host galaxy during quiescent periods, then a radio transient associated with 
such an object may be confused with a RSN or orphan afterglow if optical observations were made when 
the AGN was in a low state. An example is the source SDSS J124602.54 which has been proposed as a 



candidate optical orphan afterglow, and later confirmed to be a BL Lac ( Gal- Yam et al. 2002 ). In order to 



select out such sources in future surveys, monitoring in the radio and/or optical bands may be necessary. 



5. Discussion 

We have calculated the number of orphan radio afterglows expected to be detected in a flux limited 
survey. For the GRB parameters derived by Frail et al. (2001); total jet energy E ~ 10^^ erg and average 
beaming factor < /^^^ >= (6*^/2)"^ ~ 500, the number of sources expected to be present over all sky at any 
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given time with flux exceeding /^min is [see Eq. (jy)] 

Nn^2Q( % V^'^^T'- (15) 

Here n is the GRB rate at redshift z = 0, n = 0.1?i_icm^'^ is the number density of the ambient medium 
into which the blast wave expands, and is the energy density of the magnetic field behind the shock in 
units of the equipartition value, fi is normalized in Eq. ( p^ ) to the value obtained under the assumption 
that the GRB rate follows the redshift evolution of star formation rate ( Schmidt 2001 ). If the GRB rate 
does not evolve with redshift, yi is a factor of 8 larger than the value used in Eq. (|l5). Since afterglow 
observations generally imply > lO^'^cm"'^, and values ^ Icm^'^ are obtained in several cases [see 
discussion in the paragraphs following Eq. (|^)], our choice of n and values are conservative. Even for 
such a conservative choice of parameters we expect « 20 radio afterglows to be detected by an all sky survey 
with /lymin — 5 mjy. These sources are detected out to a distance 

<„=»2„0„=?(i|)"'(/=^)"'"Mpo, (10) 

We conducted a search for point-like radio transients with flux densities larger than 6 mJy using the 
FIRST and NVSS surveys, and discovered 25 radio afterglow candidates, listed in Table 1. The two main 
types of sources that may produce radio transients, which may be confused in our analysis with a radio 
afterglow, are radio supernovae (RSNe) and radio loud AGNs. While we have shown that our candidate 
sources are unlikely to be RSNe (see sec. §|4.2|), the possibility that most of them are radio loud AGNs 



cannot be ruled out (see sec. §4.3). Thus, our detected sources allow to set an upper Hmit to the number of 



orphan radio afterglows. Correcting for various completeness factors (see sec. §4.1), we obtained an upper 
limit of 447 all sky radio afterglows above 6 mJy, and 154 above 10 mJy. 

We have shown that for a given isotropic equivalent burst energy, E-^so — fh^E, the number of 
afterglows detected in a flux limited survey is smaller for larger beaming factor /^^^, i.e. for smaller jet 
opening angle 9 [see Eq. (^]. The upper limit we derived on the number of radio afterglow sources therefore 
implies a lower limit for the beaming factor, 

/ ■ \ 6/5 / c \ 27/20 



0.5Gpc-V"V VO-03 

Here we have used Eq. (^ with < eI^J^ >6/i1— 5 x 10^'^ erg, based on GRBs with known redshifts (Frail 
et al. 2001). This result is consistent with the value < /^^^ >= 500 inferred by Frail et al. (2001). Better 
determination of the distribution of n and ^s, based on afterglow observations of identifled GRBs, will 
improve the constraint on beaming. Alternatively, if the beaming factor is determined independently using 
other methods, then eq. ( p7| ) can be used to constrain the afterglow parameters. For instance, invoking 
< f^^ >= 500 implies 

■ * \ 27/19 / . X -24/19 



n < 5 = cm~3. (18) 

,0.03/ V0.5Gpc-V-V 

In reality the afterglow parameters in a sample of radio GRBs should have a spread, and so the letter 
condition may be taken as a constraint on the mean. 

Additional observational efforts aimed at identifying the candidates in table 1 are strongly motivated 
by two arguments. First, some of the candidates may turn out to be orphan afterglows. The fact that most 
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candidates in table 1 have no optical counterparts brighter than the POSS-II limit, implies that if some of 
them are indeed GRB radio orphans then their host galaxies are not as bright as would have been expected 
under the assumption that GRBs are associated with star-forming regions, given the distances estimated in 
Eq. (|l^) for the conservative values of n and £^b- Nevertheless, the observed number of orphan candidates 
is consistent with larger values of these parameters, as shown in Eq. (^8|) , for which the distances are larger 
and the absence of optical counterparts is not surprising. Second, further observations are also motivated 
by the fact that rejection of sources from our list in Table 1 as possible candidates would impose a more 
stringent constraint for the beaming factor. Rejection of all the sources in Table 1 would imply an all sky 
upper limit of 49 sources, and a corresponding beaming factor /j^^ > 100. 

Future radio surveys may go deeper than the FIRST and would be of comparable or better quality. The 
Allen Telescope Array (formerly the " 1-hectar radio telescope" ) will allow a full coverage of the northern sky 
every night with ~ 1 mJy sensitivity]^, and the square-kilometer-array may provide even better sensitivity^. 



A search for point-like radio transients similar to that outlined in §4.1, would probably yield a large number 
of candidates. The analysis presented in ^ predicts that the number of all sky radio GRBs above 1 mJy is 
likely to lie in the range between a few hundreds and a few thousands, depending on afterglow parameters 
[see Eq. (^5|)], provided the average beaming factor of GRBs is not significantly larger than currently 
estimated from other considerations. Since the detectable orphan population is dominated by ~ 1 yr old 
GRB remnants, which have just undergone the transition to sub-relativistic spherical expansion, these 
sources should be transient, exhibiting ~ 1/i flux decline on time scale of months. Expansion at the speed 
of light over ~ 1 yr implies that the angular size of the nearby remnants, lying at a distance ^ 100 Mpc, is 
w 1 mas. The nearby remnants may therefore be resolved using the VLBA and, moreover, since they have 
just undergone the transition to spherical expansion, their structure may show significant anisotropy (see 
Paczyhski 2001). 

RSNe and radio loud AGNs would probably contribute a large number of transients, and may render 
the identification of radio afterglows difficult. As shown in §^]^, RSNe and RGRBs can, in principle, be 
distinguished statistically by their redshift distributions. At the mJy level, the average distance of the radio 
afterglows is expected to lie between 0.1 to 1 Gpc, depending on afterglow parameters [see Eq. (|l6|)], so 
it should be possible to detect the host galaxies and obtain their redshifts. RSNe should typically lie at 
significantly shorter distances (see § |4.2| ). Detection of such bimodality in redshift distribution will greatly 
assist in identifying radio GRBs. Furthermore, the relativistic expansion of the GRB remnants also implies, 
as explained above, that unlike RSNe they may be resolved with the VLBA, possibly revealing anisotropic 
structure. Efficient rejection of flaring blazars would require follow-up observations in the radio, optical and 
X-ray bands. In particular, the orphan radio afterglows should exhibit a decay of the radio flux. 

Finally, it is conceivable that GRBs emit substantial amounts of energy in the form of gravitational 



waves and high-energy neutrinos ( 


Waxman & BahcaU 1997, 


Rachen & Meszaros 1998, 


Waxman & Bahcall 


1999, 


Waxman & Bahcall 2000 




Meszaros, & Waxman 2001 


; see 


Waxman 2001 


for review). Association 



of orphan radio afterglows with potential signals in gravitational wave and neutrino detectors may help 



increasing the signal-to-noise ratio. In a model proposed recently (Van Putten 2001, Van Putten & Lcvinson 



2002| ) , the GRB results from extraction of the rotational energy of a Kerr black hole through its interaction 



with a magnetized torus. This model predicts that the major fraction of the rotational energy available will 
be emitted in the form of gravitational waves in the band 0.5 - 1.5 KHz, for a reasonable range of black 



5 = 



see tttp:/ /astron. berkclcy.edu/ral/ 



^http:/ /www. nfra.nl/skai/ 
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hole masses. It further proposes that measuring the product of the total energy carried by the gravitational 
waves and the corresponding frequency, which reflects the compactness of the system, can provide an 
existence test for Kerr black holes. However, such a measurement requires a knowledge of the distance to 
the source. Since the gravitational wave emission is isotropic, the chance of associating a LIGO source with 
a GRB event is small if GRBs have a beaming factor as large as currently believed. Nevertheless, if this 
model is correct, then LIGO sources should be associated with the orphan radio afterglows that should 
appear several months after the burst of gravitational radiation. At the distance within which LIGO can 
detect the predicted signal (~ 100 Mpc), the radio source is expected to be very bright (see eq. jl^) and 
can be monitored for at least several years. Even though the angular resolution of the gravitational wave 
telescopes is expected to be quite poor, the chance that such a bright transient would coincide with a 
LIGO source should still be very small, and so such an association would be statistically significant. Our 
analysis for instance yielded an upper limit of about 10^^ transients per square degree above 10 mJy. The 
association of a radio afterglow with a LIGO source would also provide a means to identify the host galaxy, 
given the superior resolution of the radio telescopes, and hence a distance measurement for the LIGO 
source. 

We thank K. Weiler, A. Loeb B. Paczyhski and R.D. Blandford for useful discussions. EW is partially 
supported by AEG grant and a Minerva grant, and is an incumbent of the Beracha Career development 
Ghair. 
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Table 1: A list of radio- transient candidates. 



H.A. 


Dec. 


I'lKbl radio 


Quality 


Optical counterparts and comments 


(J2000) 




ilux [mjyj 






12:56:32.34 


-05:47:39.7 


9.8 


2 


No optical source 


13:07:13.39 


f c\n r\r\ a 

-05:27:09.4 


6.6 


2 


No optical source 


12:26:07.68 


-1-02:02:52.4 


12.9 


1 


No optical source 


12:28:47.38 


+02:03:12.8 


34.5 


3 


No optical source 


12:29:02.49 


+02:31:03.1 


9.6 


2 


No optical source 


12:26:50.43 


+02:39:32.6 


6.8 


1 


Blue point source with it ~ 19 mag 


10:08:02.94 


+07:33:18.2 


6.5 


3 


No optical source 


12:30:23.41 


+ 11:31:39.0 


7.6 


2 


Point source with ft ^ 18 mag 


LZ.Zo.oZ.bZ 


+ iz:z:o:UU.o 


I .1 


A 


Nucleus of an i? ~ 16 galaxy 


12:29:00.02 


+12:48:18.8 


6.2 


2 


No optical source 


12:15:50.24 


+ 13:06:54.0 


9.7 


2 


Located in an NGC4216 galactic arm 


08:21:50.18 


+17:46:16.3 


6.0 


2 


No optical source 


11:43:55.34 


+22:10:20.4 


6.2 


2 


No optical source 


16:52:03.08 


+26:51:39.9 


6.3 


2 


Coincident with the PSR J1652+2651 


13:31:14.55 


+30:25:57.8 


10.6 


3 


No optical source ^ 


13:31:18.21 


+30:26:01.1 


10.8 


3 


No optical source ^ 


13:31:01.86 


+30:33:00.8 


13.1 


3 


No optical source ^ 


13:31:14.42 


+30:55:39.2 


8.7 


2 


No optical source ^ 


13:31:17.11 


+31:01:53.2 


7.4 


2 


No optical source ^ 


13:24:20.39 


+31:41:37.2 


9.2 


1 


0.2' from an i? ~ 16 galaxy 


13:25:17.11 


+31:51:07.6 


9.4 


2 


No optical source 


17:20:59.90 


+38:52:26.6 


8.2 


2 


No optical source 


15:22:48.69 


+54:26:44.1 


6.6 


2 


Point source with i? ~ 19 


10:48:48.92 


+55:15:08.7 


6.1 


2 


Blue point source with B ^ 21 


08:36:00.07 


+55:55:17.9 


10.7 


2 


No optical source 


Notes: 



" radio images quality: 1 - probable radio transient, 2 - suspect detection, 
3 - probably an artifact 

^ group of 'sources' near bright radio source 3C286 - probably artifacts. 



